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The connection between statics and dynamics of spin glasses 
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We present results of numerical simulations on a one-dimensional Ising spin glass with long-range 
interactions. Parameters of the model are chosen such that it is a proxy for a short-range spin glass 
above the upper critical dimension (i.e. in the mean-field regime). The system is quenched to a 
temperature well below the transition temperature Tc and the growth of correlations is observed. 
The spatial decay of the correlations at distances less than the dynamic correlation length ^{t) agrees 
quantitatively with the predictions of a static theory, the “metastate”, evaluated according to the 
replica symmetry breaking (RSB) theory. We also compute the dynamic exponent z{T) defined by 
^{t) oc t^C(T) that it is compatible with the mean-field value of the critical dynamical 

exponent for short range spin glasses. 


Experimental measurements on a system at finite- 
temperature involve a time average. Dynamics is harder 
to calculate than statics, so, in theoretical work, the time 
average is usually replaced by a static calculation using 
statistical mechanics in which one sums over all con¬ 
figurations with the Boltzmann probability distribution. 
Most systems are ergodic, so theory agrees with experi¬ 
ment even though different averages are performed. One 
situation where more care is needed is that of a phase 
transition where symmetry is spontaneously broken. A 
simple case is the Ising ferromagnet, which has just two 
ordered states below the transition temperature Tc, the 
“up” spin state with a net positive magnetization, and 
the “down” state. On cooling the system will choose one 
of these symmetry-related states and acquire a non-zero 
magnetization. However, the Boltzmann sum will (un- 
physically) include both the up and down states and give 
zero net magnetization. The “up” and “down” states 
in the Ising ferromagnet are called “pure states” in the 
literature [l(]. In pure states, correlation functions have 
a “clustering” property which means that “connected” 
correlation vanish at infinity, i.e. 

hin i{S.S,)-{S.){S,)) = 0. (1) 
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This does not occur in the Boltzmann average since the 
second term is zero, so the combination of “up” and 
“down” is not a pure state; rather it is a “mixed state”. 

However, in spin glasses, which have disorder and 
“frustration”, the situation is much more complicated. 
Dynamically, below the spin glass transition temperature 
Tc a macroscopic system is not in thermal equilibrium be¬ 
cause relaxation times are much too long. Rather, in a 
typical experiment the system is quenched from a high 
temperature to a temperature below Tc and the subse¬ 
quent dynamic evolution of the system is observed. For 
statics, the state (or states) of thermal equilibrium are 
very complicated and are not related to any symmetry. 
As for the ferromagnet we would like to find a static cal¬ 
culation which will predict the experimental behavior, at 
least to some extent. In this paper we show quantitatively 
that the theoretical idea called the “metastate” 
combined with the technique of “replica symmetry break¬ 
ing” (RSB) an provides such a description for spin 


glasses, at least in high dimensions, where the critical 
behavior is described by mean field theory. 

Since the clustering property in Eq. © is convenient 
one would like to also describe spin glasses in terms of 
pure states. This can be done (in principle) by taking 
a very large system, applying some boundary condition 
on it, and studying the correlations in a relatively small 
window far from the boundary P, Q . This procedure is 
repeated for many different boundary conditions. 

The question of whether there are many pure states or 
just one (a time-reversed pair in the absence of a mag¬ 
netic fiel^ in spin glasses has been very controversial 
i a a a, If there are many, one needs to do some 
sort of statistical average over them, which is called a 
“metastate”, for which different but equivalent formula¬ 
tions have been given by Newman and Stein 0 and by 
Aizenman and Wehr (AW) Q ■ In Hi® metastate, one 
considers a scale M, intermediate between the window 
size W and the system size L. The metastate-averaged 
state (MAS) is obtained by computing correlation func¬ 
tions in the window in which an average is performed 
not only over the spins but also over the bonds in the 
“exterior” region between M and L [1, H, |^ . 

For the infinite-range model of Sherrington and Kirk¬ 
patrick (SK) 0 , Parisi’s exact solution [3, solved by 
RSB, has many pure states in a sense that was clarified 
later by Newman and Stein 0, see also Read 

The critical behavior of a realistic ^n glass is expected 
to be the same as that of the SK model in dimen¬ 
sion d greater than the “upper critical dimension”, du , 
which is equal to six. However, this does not necessarily 
mean that the RSB description of the g)in glass phase 
below Tc also applies for d > 6 [1, 0, 111- Nonetheless 
Read [1] has computed the spatial fluctuations in a finite¬ 
dimensional model below Tc, assuming mean-field (Gaus¬ 
sian) fluctuation^ and the metastate description coming 
from Parisi’s [3, Q RSB solution of the SK model. Spin 
correlations are found iiniiii to decay with a power of 
the distance, due to the averaging over many pure states 
(which are unrelated by symmetry) in the metastate, i.e. 

{SiSj)MAs « 1/g/ with as=d-A, (2) 

where “s” refers to “static, “MAS” stands for metastate- 
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averaged state, and sites i and j are in the window far 
from the boundary. The result in Eq. @ had been ob¬ 
tained earlier in Ref. [H from an RSB calculation work¬ 
ing in the zero overlap sector. 

We emphasize that the calculation leading to Eq. © is 
a static one. Is it possible to relate it to experiments (or 
numerical simulations), which concern (non-equilibrium) 
dynamics! Many simulations [§, [l3l - [l^ have been car¬ 
ried out in which a spin glass is quenched to below 
and the resulting dynamics analyzed. It is found that 
fluctuations reach a steady state on length scales smaller 
than a dynamic correlation length ^(t) which is found, 
empirically, to grow with a power of t like 

m oc , (3) 

where the non-equilibrium dynamic exponent z{T) is 
found to vary, roughly, like 1 /T and becomes close to the 
critical dynamical exponent, Zc, for T = Tc, [l3l - [l^ 

llz{T)^{TlT,)Zc. (4) 

At distances less than ^(t) correlations are observed to 
fall off with a power of distance leading to the following 
scaling prediction 

c,{n,,t) ^ [(5.(t)^,(t))2] = ^ / , (5) 

where “d” refers to “dynamic”. Here the thermal average 
squared, is performed by simulating two copies of 

the system with the same interactions, initialized with 
different random spin configurations. Spin configurations 
of the two copies at the same time are used in Eq. (O. 
Use of two copies provides an unbiased estimate of this 
thermal average. The second average, [•••], is over the 
bonds. We will also average over all pairs of sites a given 
distance r apart. Eor <C C{t) one finds f{x 0) = 
const, so 

C 4 {rij,t) oc l/r^f {rij < ^{t)). (6) 

Clearly, the non-equilibrium dynamics is generating a 
sampling of the pure states. To our knowledge. White 
and Eisher were the first to point out the similar¬ 
ity of this sampling and the metastate average for stat¬ 
ics. They use the term “maturation metastate” to de¬ 
scribe the ensemble of states generated dynamically on 
scales less than ^{t) following a quench, and “equilib¬ 
rium metastate” for the static metastate discussed ear¬ 
lier. Here we will use terms “dynamic” and “static” to 
describe these two metastates. Subsequently Manssen 
et al. [2l| emphasized the similarity between the two 
metastates and suggested that they might actually be 
equivalent, in which case as in Eq. m would equal ad 
in Eq. ®. The rationale behind this hypothesis is that 
thermal fluctuations of the spins outside the window at 
a distance ^{t) and greater, which are not equilibrated 
with respect to spins in the window, effectively gener¬ 
ate a random noise to the spins in the window which is 


similar to the random perturbation coming from chang¬ 
ing the bonds in the outer region according to the AW 
metastate. It would also be very useful if the metastates 
were equivalent because then a theory of the statics of 
spin glasses would give results corresponding to experi¬ 
ments, which are a time average, as is the case for simpler 
systems with a phase transition like ferroma gne ts. 

For the three-dimensional spin glass. Refs. [2^ and 
have shown that an equilibrium calculation in the zero 
spin-overlap sector gives a power-law decay for the spin 
correlations, as in Eq. with a value of as consistent 
with that obtained from dynamics following a quench in 
Ref. 0. These are both numerical results. Here we wish 
to consider the mean field limit, d > 6, because there is 
an exact analytic result Ofs = d — 4, in RSB theory iEl 
to compare with. Unfortunately it is difficult to carry out 
useful Monte Carlo simulations for a spin glass below 
in more than six dimensions, because the number of sites 
in a region of linear size increases so fast, as that the 
range of £ that can be studied is very limited. The only 
calculation of the exponent ad in the mean field region 
that we are aware of is that of Ref. who studied d = 6 

but only at T = Tc and so these results are for the critical 
point rather than the spin glass phase. 

However, it has been pointed out that a class of models 
in one-dimension, with long-range interactions falling off 
with a power of distance can serve as a proxy for short- 
range models [2B - IT7I| in a range of dimensions including 
high dimension. The Hamiltonian is 

n = -J2j,,s.s,, (7) 

where the sites i = 1, 2, • • • , A^ he on a one-dimensional 
chain with periodic boundary conditions, the Ising spins 
Si take values ±1, and the interactions Jij are indepen¬ 
dent random variables with whose distribution has mean 
and variance given by 

[J.,] = 0, [4] oc 1/4/, (8) 

in which tr is a parameter which can be varied. To incor¬ 
porate periodic boundary conditions it is convenient to 
put the sites on a ring and define Rij to be the chord dis¬ 
tance between i and j, i.e. Rij = {N/tt) sin(7r|i — j|/A^), 
whereas the distance along the ring is = \i — j\ if 
I* — j\ < N/2 and = N — \i — j\ otherwise. The 
bonds are generated, and the constant of proportional^ 
in Eq. ([8]) fixed, in the following way due to Ref. (^ . 
We choose an average coordination number Zb, which we 
take here to be Zb = 6. We choose a site i at random and 
then a site j with a probability C/Rf/j where C is the 
normalization constant. If there is already an interaction 
between i and j repeat until a pair {i,j) is found which 
has not occurred before. Then assign an interaction be¬ 
tween i and j chosen from a Gaussian distribution with 
mean zero and standard deviation unity. Repeat this 
Nzb/2 times, so there are Nzb/2 interactions. 

Varying tr is argued to be analogous to changing the 
dimension d of a short-range model |26l | . In the mean field 
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FIG. 1: Data for the correlation function C 4 (r, t), defined in 
Eq. (O, as a function of r = |i —for a range of sizes between 
N — and 2^®. The data is averaged over about 1000 bond 
configurations. It is also averaged over times between 2^^ and 
half that value. There are clearly strong finite-size effects but 
the data seems to have converged for the largest sizes, at least 
up to distances of order 10"^. 


regime (for the short-range case, d > c?„ = 6), a precise 
connection can be given between a and an equivalent d, 
namely 

d = 2/(2a-l), (9) 

and so, for the long-range model, the mean field regime 
is 1/2 < a < 2/3. 

The connection between critical exponents of the 
short-range and corresponding long-range models has 
been discussed systematically in Ref. [^, where it is 
noted that an exponent of the short-range model in d 
dimensions is d times the corresponding exponent of the 
equivalent one-dimensional long-range model. Hence, to 
get the exponent = d — 4 in the static metastate for 
the long-range model we divide by d and, since we work 
in the mean held regime, use Eq. © to relate d to a. 
This gives 

as = 3 — 4a (long-range model). (10) 

Here we focus on one value in the mean-held regime, 
a = 5/8, which corresponds to d = 8 according to Eq. ([HI). 
Using standard hnite-size scaling analysis we hnd that 
Tc = 1.85(2) for this model with zt = 6. We need to 
work well below Tc so that our data is characteristic of 
the ordered phase and does not also incorporate critical 
huctuations. We take T = 0.4Tc = 0.74. We have pre¬ 
ferred to focus the numerical effort, which is substantial, 
on one temperature in order to get the best quality data 
for the largest possible range of sizes. 



FIG. 2: Data for the correlation function for the largest 
size N — 2^® as a function of r for different times between 
1 — 2"^ and 2^'^. A gradual crossover can be seen between two 
power laws. At large t and small r, C 4 {r,t) oc with 

a?d = 3 — 4 (t (which is also the value of as expected from 
the static metastate using RSB theory). This is the region 
in which the data has reached a steady state. At large r and 
small t one finds C 4 {r,t) oc l/r^*^ which is just the average of 
the square of the interactions Jij. 


We quench the system from infinite temperature to 
T = 0.74 at time t = 0 and follow the evolution of the 
system using Monte Carlo simulations. We measure spin 
correlations, averaging them for times between 2^ and 
2^“^, for integer A: up to a maximum value. For the 
largest sizes this was k = 14. We hnd that hnite-size 
effects are very large and we need to study enormously 
large sizes. We therefore take a range of sizes which also 
increases geometrically, = 2^ up to £ = 26. We also 
average over about 1000 samples (the precise number de¬ 
pending on size). 

Figure [T] shows our data for for the correlation function 
C 4 (r, t), dehned in Eq. ([5]), as a function of r = \i — j\ at 
t = 2^^ for different sizes. Despite the strong hnite-size 
effects the data seems to have converged for the largest 
sizes at least for the range of distance presented. 

Having established that the largest size, N = 2^®, is 
large enough to eliminate hnite-size effects for the range 
of r and t being considered we now discuss the data for 
this size in more detail. Figure [2] shows data for C' 4 (r, t) 
at different times as a function of r. It is expected to have 
the scaling form shown in Eq. ([5]). For short range mod¬ 
els the scaling function f{x) decays exponentially at large 
X because the correlation function falls off very rapidly 
once r is greater than the dynamic correlation length. 
However, in the present model we have interactions of ar¬ 
bitrarily long range which give a “direct” contribution to 
the correlation function at large distances. Since Ci(r,t) 
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FIG. 3: The main part of the figure is a scaling plot of our 
data for the largest size N — 2^® at T = 0.74 according to 
Eqs. ([5} and (|3}. We assume a pure power law for ^{t), namely 
with ^{t) = . The data collapses well with a dynamic 

exponent 2 (r) = 1.4. The inset is a plot of the effective 
exponent a^ff, the slope of the curves in Fig. [2] obtained by 
differentiating a spline fit. The lines are quadratic fits to the 
data for intermediate r (7 < r < 255). One sees corrections 
to the parabolic hts at very small distances, r < 3. 


involves the square of the spin-spin correlation function, 
and is averaged over the interactions, the direct contribu¬ 
tion should be proportional to [Jjy], which, according to 
Eq. dSI), is proportional to oc (= for a = 5/8). 

Fig. m follows this behavior at short times and large dis¬ 
tances, see the dotted line. 

By contrast, at small r and large t, where r f (t), the 
data for different times collapses and is consistent with 
a decay (= for a = 5/8), see the dashed 

line in Fig. [21 To estimate better the slope at large t 
and small r we plot in the inset to Fig. O the effective 
exponent Oeff, the slope of the data in Fig. [21 as a function 
of r for different times. The curves are quadratic fits for 
intermediate r {7 < r < 255). The intercepts of the fits 
approach —0.5 for r —!■ 0 at large t. Hence, according to 
Eq. (ini), we have = 3 —4 ct (or at least very close to it.) 
However, this is precisely equal to as, the corresponding 


exponent from the static metastate according to RSB 
theory as shown in Eq. m- Hence we see that, in the 
mean field regime, the static and dynamic metastates 
appear to agree and the description appears to be that 
of RSB. The latter agrees with some other studies [s^, 
and is also implied by those, such as Refs. [H, [2^, which 
argue that RSB holds even below six dimensions. 

The main part of Fig. [3] shows a scaling plot of our 
data for the largest size according to Eqs. m and 
The data scales well and indicates 2;(0.4Tc) = 1.4(2). For 
short-range models z{T) is found to obey Eq. (jj]). If we 
assume the same here then z{Tc) = Zc = 0.56(8). To 
translate this value for Zc, the critical dynamical expo¬ 
nent, into the exponent for the equivalent short-range 
model, we multiply by d (= 8), as discussed above, so 
our estimate for the critical dynamical exponent of the 
d = 8 short-range spin glass is = 4.5(6) (d = 8). This 
model is in the mean held region (d > 6) for which the 
dynamical exponent is found to be Zc = 4 MM- Our 
result is consistent with this. 

To conclude, we have shown quantitatively that the 
non-equilibrium dynamics following a quench of a model 
which is a proxy for a short-range spin glass in dimen¬ 
sion d > 6 is given, in the steady-state regime where 
the distance is less than the non-equilibrium correlation 
length, by the analytic result for the static metastate cal¬ 
culated according to RSB theory. This suggests that (i) 
RSB theory applies to spin glasses above the upper crit¬ 
ical dimension, d„ = 6, and (ii) the dynamic and static 
metastates are equivalent (at least in this region). If the 
latter is true in general it would provide a great simpli- 
hcation in the study of spin glasses. 
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